Radio frequency pulsed-gate charge spectroscopy on coupled quantum dots 
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Time-resolved electron dynamics in coupled quantum dots is directly observed by a pulsed-gate 
technique. While individual gate voltages are modulated with periodic pulse trains, average charge 
occupations are measured with a nearby quantum point contact as detector. A key component of 
our setup is a sample holder optimized for broadband radio frequency applications. Our setup can 
detect displacements of single electrons on time scales well below a nanosecond. Tunneling rates 
through individual barriers and relaxation times are obtained by using a rate equation model. We 
demonstrate the full characterization of a tunable double quantum dot using this technique, which 
could also be used for coherent charge qubit control. 

PACS numbers: 03.67.-a, 06.60.Jn, 73.63.Kv, 84.40.-x 



I. INTRODUCTION 

Solid state quantum information processing based on 
single electrons in coupled quantum dots (QDs) demands 
very precise control of the quantum states involved (TJ |2] . 
On the one hand, it is necessary to perform many qubit 
operations within the coherence time. On the other 
hand, to address quickly specific qubits in a large ar- 
ray and detect their response poses a serious experimen- 
tal challenge. On the detection side, single-shot readout 
of charge states of QDs has been demonstrated up to a 
bandwidth of 10 MHz which is limited by the detector 
noise [SHTT]. Shorter time scales can be resolved using 
a time-averaged detection while an excitation pulse cy- 
cle is continuously repeated |T2HT6|. In this way, charge 
coherence times below 1 ns have already been resolved 
in transport spectroscopy measurements [151 US] where 
the chemical potential of one lead of a double QD was 
pulsed while the average current through the double QD 
was measured. However, this technique, based on pulsing 
ohmic contacts, will fail in more complex devices which 
consist of a larger number of coupled QDs with only 
few ohmic contacts (a more realistic circuit for quantum 
information processing). The accessibility of individual 
QDs, needed for the initialization or readout of specific 
qubits, requires more involved techniques such as high- 
bandwidth locally resolved charge spectroscopy. In this 
scenario, several close-by gates of a complex nanostruc- 
ture can be individually pulsed while the conductance of 
a capacitively coupled quantum point contact (QPC) [17] 
or alternatively a detector QD is measured. A main ob- 
stacle for high-bandwidth charge spectroscopy is cross- 
talk between adjacent gates, that is, leakage of a volt- 
age pulse originally applied to one specific local gate into 
other components of the nanostructure. 

We developed a high-bandwidth pulsed-gate technique 
with square pulses much shorter than a nanosecond. 
Here, we use it to study time-resolved energy relaxation 
and quantum-mechanical tunneling of a single electron in 
a double QD. To avoid the disadvantages of current spec- 



troscopy, we perform charge spectroscopy using a nearby 
QPC as a detector. A similar pulsed-gate technique has 
already been used for the control of spin qubits [T8ti20) . 
However, compared to our measurements, the pulse du- 
rations and repetition periods reported there were much 
longer, hence, imposing smaller challenges to radio fre- 
quency (rf) components. In addition, these earlier and 
following works relied on a qubit initialization via charge 
exchange with the leads, which involves pulsing between 
three different charge configurations of the double QD. In 
contrast, we pulse between only two configurations and 
work at a constant overall charge. In this case, charge 
exchange with the leads is no longer required, and the 
tunnel coupling to the leads can be tuned to be very 
weak. We expect this effective decoupling from the bath 
of two-dimensional electrons to ultimately lead to longer 
coherence times. Here, we demonstrate a characteriza- 
tion of the double QD parameters relevant for the control 
of a charge qubit. We use the same pulsed-gate technique 
which will be employed to actually operate qubits, hop- 
ing to increase the overall efficiency of a future quantum 
computer. Introducing a rate equation model allows us to 
obtain from our data tunneling rates of individual barri- 
ers and energy relaxation times. Moreover, we are able to 
distinguish between possible charge relaxation channels, 
which can be tuned via gate voltages. This is especially 
important for future qubit applications where specific re- 
laxation channels might be temporarily opened for ini- 
tialization or readout, while long energy relaxation times 
are required during coherent manipulation. 



II. NANODEVICE AND SETUP 

Mapping the applied pulses with a short rise time 
(~ 70 ps) to a time-dependent electrical potential at the 
double QD requires high rf bandwidth at cryogenic tem- 
peratures. Our gate geometry and rf sample holder are 
designed to minimize the detrimental cross-talk and op- 
timize impedance matching. At the same time, using 
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FIG. 1: (a) A scanning electron micrograph of a device 
nominally identical to the investigated one. Details in the 
main text, (b) The radio frequency sample holder with four 
impedance-matched micro strip lines and 22 low bandwidth 
connections. The central gold square is electrically shorted 
to the copper plated backside of the circuit board. The sam- 
ple would be mounted on the central gold square and wire 
bonded. 

rf wiring in our dilution refrigerator, we have extended 
the bandwidth of our measurement to beyond 10 GHz 
(presently limited by the use of micro miniature coax 
connectors). 

Figure [T^ shows a scanning electron micrograph of 
the surface of the double QD device based on a 
GaAs/AlGaAs heterostructure. It contains a two- 
dimensional electron system at a temperature of Tq\ c:^ 
130 mK. The metallic gates can define electrostatically a 
few-electron triple QD [21, 22 . However, in the present 
experiment, we only define a double QD (labeled B and 
C) and implement one of three QPCs as a charge detec- 
tor [17]. Therefore the three darkened gates (a, 61 and 
^^1) in Fig. ^ are connected to the electrical ground. 
All ohmic contacts (labeled I - VI) are also grounded 
with the exception of contact IV, which is biased with 
Vqpc — — 50/iV. The resulting dc current, /qpc, flows 
through the QPC to the grounded contact III , while cur- 
rent flow through the double QD is negligibly small even 
under resonance conditions, /qpc depends on the charge 
state of the double QD. We reduce back-action of this 
charge sensor onto the double QD by limiting the detec- 
tor sensitivity (at still high contrast) [23 . In detail, we 
only apply a low voltage Vqpc and operate the QPC at 
a small conductance of Gqpc ^ 0.5 x 2e^ /h. 

Gates tl, /3, t2 and 7 are connected to rf coaxial cables. 
The dc voltage applied to these gates can be modulated 
via high-bandwidth bias tees (Agilent 11612 A) operated 
at room temperature. Pulse trains are produced with a 
pattern generator (Agilent 8 1134 A) that provides a min- 
imal rise time of 70 ps. A main technological challenge is 
the precise transmittance of the pulse shapes generated 
at room temperature to individual gates of the nanoscale 
device at low temperature. When choosing rf cabling, 
a compromise has to be made between a small thermal 
conductivity and a low rf damping. We use stainless 
steel coaxial cables with a silver-plated inner conduc- 




FIG. 2: Charge stability diagrams of the double QD. Plotted 
in gray scale as a function of gate voltages and (see 
Fig.[ljL) is the transconductance d/qpc/dl/g, which has been 
obtained by numerical differentiation of the dc current /qpc- 
d/Qpc/dV/3 :^ at the gray background, d/gpc/dy^ < on 
darker lines and d/Qpc/dV/g > at brighter regions. Pairs 
of numerals or 1 depict the number of electrons charging 
the QDs B,C in the ground-state configuration, (a) No pulses 
applied, (b) is modulated with a symmetric continuous 
square wave signal, for which the waiting time tw between 
two pulses equals the pulse duration tp, at a period of /q"^ = 
tw + = 554 ps and an amplitude of = 5 mV (see inset 
and main text). 



tor from room temperature to T c:^ 4.2 K and supercon- 
ducting coaxial cables at lower temperatures. For ther- 
mal anchoring of the inner conductor of the coaxial lines 
at various temperatures, we have developed impedance- 
matched gold strip lines on sapphire sheets with insertion 
losses below 0.4 dB. These devices are virtually resonance 
free up to 12 GHz. Our broadband sample holder ^4] 
(Fig. ^p) contains four micro strip lines which are op- 
timized - with the help of numerical calculations using 
the commercial software SONNET - for transmission of 
short pulses and minimal cross-talk [8 . The sample is 
mounted on the central gold square and electrically con- 
nected via short bond wires. 



III. STABILITY DIAGRAM AND PULSE 
CALIBRATION 

Figure [2^ is a plot of a section of a charge stability 
diagram of our double QD, where it contains between 
zero and two electrons, as a function of gate voltages 
and Vp. The transconductance dloiPc/dV^ (gray scale), 
obtained by numerical differentiation, clearly shows four 
charging lines (dark) separating stable ground-state con- 
figurations with zero or one electron in QD B and C, 
where pairs of numerals or 1 depict the number of 
electrons charging the QDs B and C. The charge re- 
configuration line (white) which separates configurations 
1,0 and 0,1 appears as a double line because of ex- 
ternal noise or back-action from the biased QPC [23] . 
The measurement is repeated in Fig. while gate P 
(y axis) is modulated with a continuous square pulse 
train at a repetition frequency of /o = 1.805 GHz and 
an amplitude oi = 5mV. In this measurement, the 
pulse duration, tp, equals the waiting time between two 
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pulses = = 277 ps (see inset of Fig. where 
+ = fo^' We call this a symmetric pulse train in 
contrast to an asymmetric pulse train with 7^ ^p- 

As a result of the symmetric pulse train, the stability 
diagram in Fig. (2)3 splits into two copies shifted along the 
y axis by the effective pulse amplitude Vp 5 mV. This 
steady- state result is predicted (for tp = t^) by ther- 
modynamics, even though our repetition rate is much 
shorter than the relevant relaxation rates - in contrast 
to earlier pulsed-gate measurements [18 . Using a cal- 
ibration via nonlinear transport measurements, Vp can 
be converted to an energy shift of the electronic spec- 
trum of the QDs which is AE 625 /ieV. This value is 
large compared to the energy scale hfo ^7.5 /ieV cor- 
responding to the pulse train repetition frequency. For 
the opposite case of hfo > AE, effects like coherent de- 
struction of tunneling would be expected [25]. For the 
limit of hfo ^ photon-assisted tunneling [2F would 
influence the outcome of similar measurements. 

Comparison of the effective pulse amplitude Vp 
5mV with the actually applied modulation amplitude 
AV = 6mV quantifies the losses during the transmis- 
sion of the pulses from the instrument at room temper- 
ature to the gate. The corresponding power attenuation 
of 201og(Vp/AV)dB -1.6 dB can be traced back to 
the damping of the coaxial cables. 



IV. RADIO FREQUENCY SAMPLE HOLDER 

To elucidate the importance of an rf sample holder, in 
this section we compare its rf properties with a nonop- 
timized sample holder combined with a standard chip 
carrier. Pulsed-gate spectroscopy with the nonoptimized 
setup suffers from reflections and strong cross-talk be- 
tween adjacent leads. This is illustrated in Fig. [3^ 
which compares network analyzer measurements of a 
nonoptimized sample holder (blue dashed line) and a rf 
sample holder (black solid line) engineered for optimal 
impedance matching and minimal cross-talk. Plotted is 
the logarithm of the power, reflected from a single open- 
ended lead (for the rf sample holder shown in Fig. [l]3, 
on of the four strip lines at the top), in the frequency 
range < / < 12 GHz. A perfect open-ended lead would 
show full reflection (0 dB), which is almost reached in 
a numerical calculation with SONNET for our rf sam- 
ple holder (red filled circles). Our rf setup is not perfect 
but its reflection losses are already much smaller than for 
the nonoptimized device. The strong oscillations of the 
reflection coefficient are caused by impedance mismatch 
resulting in resonances, phase shifts and radiation losses. 
The inset of Fig. ^ shows the calculated power spec- 
trum of a typical pulse train (tp = = 277 ps, rise time 
70 ps). It suggests that most power is contained in the 
spectral range < / < 5 GHz, where the measured reflec- 
tion properties of our rf sample holder (solid black line in 
Fig. |3^) are almost perfect. The imperfections at higher 
frequencies can nevertheless influence the pulse shape. 




FIG. 3: (a) Reflection coefficient (Sll in dB) of a single open- 
ended lead of an rf sample holder (black solid line) , optimized 
for minimal cross-talk and impedance matching, compared to 
that of a nonoptimized sample holder (blue dashed line) as a 
function of frequency. The (red) closed circles are a numerical 
simulation of the rf sample holder using the commercial soft- 
ware SONNET, (b) Transmission (S12 in dB) through the two 
sample holders as a function of frequency. The transmission is 
measured between two adjacent but electrically isolated gates. 
The inset shows the calculated Fourier transform of a pulse 
train with tp = tw = 277 ps as a function of frequency, where 
P is the power in arbitrary units, (c) A section of the double 
QD stability diagram, as in Fig.|2]3, but using a nonoptimized 
sample holder while a symmetric pulse train with a repetition 
frequency of /o = 300 MHz is applied at gate 7. (d) Mea- 
surement as in c) but at a slightly larger repetition frequency 
of /o = 400 MHz (nonoptimized sample holder). The two 
frequencies are marked in b) by vertical arrows. 

For instance, strong cross-talk between the control gates 
can lead to even more serious pulse deformation. This 
cross-talk can be quantified by measuring the transmis- 
sion between adjacent but electrically not interconnected 
leads. Figure shows typical transmission curves for 
the two sample holders in the relevant frequency range. 
The measured cross-talk of the rf sample holder is fairly 
small and almost coincides with the numerical calcula- 
tion done with SONNET (red circles). The nonopti- 
mized sample holder (blue dashed line) shows a much 
larger transmitted signal — especially in the frequency 
range of the largest power density (see inset) — and, in 
addition, strong oscillations. 

Figures jsj) and |3]i demonstrate how the transmission 
influences the cross-talk in pulsed-gate experiments. The 
two plots show the same stability diagram section of 
our double QD measured with the nonoptimized sample 
holder under application of a symmetric pulse train to 
gate 7 {x axis). While the nominal pulse shape with 
AU = 4mV and a rise time of 70 ps is identical for 
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FIG. 4: Charge stability diagrams of the double QD as in Fig. 
[2] but under application of an asymmetric pulse train (tw ^ 
^p)- and are both modulated in opposite directions. 
A square pulse of duration tp = 1 ns is periodically applied 
after the waiting time of tw = 50 ns. Effective amplitude and 
direction of the pulses are indicated by white arrows. The 
pulse amplitude Vp is increased from a) to c). 
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both cases, the repetition frequencies are slightly differ- 
ent, namely /o = 300 MHz for Fig.jSf) and /o = 400 MHz 
for Fig. |3]i. As explained above (compare Fig. 
application of the pulse train results in a splitting of 
the more vertical charging line while the more horizon- 
tal one is only broadened. The broadening is, how- 
ever, larger for the slightly smaller repetition frequency 
(Fig. [sj)). The difference between the measurements 
in Figures [S]^ and |3]i can be traced back to the trans- 
mission spectrum, which displays stronger cross-talk for 
/o = 300 MHz (arrows in Fig. [sJd). Compared to mea- 
surements using an rf sample holder (Fig. |2]3), the us- 
age of the nonoptimized setup causes all charging lines 
to be broadened, which suggests imperfect pulse shape 
along with strong cross-talk. In summary, for quantita- 
tive charge spectroscopy using pulsed-gate measurements 
with subnanosecond time scales, the usage of a sample 
holder optimized for rf frequencies is essential. 



FIG. 5: (a) Qualitative explanation of the electron dynam- 
ics for a asymmetric pulse train, (a-i) Chemical potentials 
of the QDs and leads (/ii) in the middle of the black- white 
line in Figs. [4^ ~[4]^- Vertical lines symbolize tunnel barriers. 
In the ground-state configuration 0,1, one electron occupies 
QD C. (a-ii) Situation while the double QD is pulsed onto 
(the middle of) the charge reconfiguration line. The chemical 
potentials of both QDs are degenerate, which allows the elec- 
tron to tunnel between the QDs with rate Fi. (a-iii) After the 
pulse, the electron occupies either QD B (iii-b) or is back in 
QD C (iii-a). (a-iv) Two possible relaxation channels (com- 
bined rate from the excited state (with QD B occupied) 
into the ground- state. One channel (iv-a) features an inter- 
dot tunneling process in combination with energy relaxation. 
The other channel (iv-b) involves two resonant QD-lead tun- 
neling processes and energy relaxation in the leads, (b) The 
sketch shows the steady-state probability of the (unpulsed) 
ground-state being unoccupied, P(t), as a function of time 
for arbitrarily chosen Fi and F2. See main text for details. 
The periodic pulse train is also indicated where tp is the pulse 
duration and tw is the waiting time between pulses. 



V. PULSE SEQUENCE 

Figures [4^ -|4j3 show measurements of the same region 
of the stability diagram as in Fig.[2]but now for asymmet- 
ric pulse trains with tp <C t^ (^p = 1 ns and t^ = 50 ns). 
The pulses are applied in the direction of the white arrows 
parallel to the charging line of QD C. For this purpose, 
the two plunger gates /3 and 7 are modulated simultane- 
ously and with opposite sign {V^ is increased while is 
decreased) [27 . From Figs. |4^ - the pulse amplitude 
is stepwise increased. For the applied asymmetric 
pulse train, a naive thermodynamic model would predict 
the average charge occupation to almost correspond to 
the unpulsed ground-state configuration — according to 
the local potential present between pulses. Within this 
model, we would expect to observe the same stability dia- 
gram as the one in Fig.|2|L (with no pulses applied). This 
is indeed the case in Figs. [4^ -[4]3, with the exception of 
an additional line (black-white) which corresponds to an 
average nonequilibrium charge configuration. The black- 
white line copies the charge reconfiguration line (white 
double line) in terms of having the same length and being 



parallel to each other. The distance between the black- 
white line and the charge reconfiguration line equals the 
effective pulse amplitude l/p. Hence, the local potential 
during the short pulses is identical to that on the charge 
reconfiguration line (in a measurement without pulses). 
Note that the black-white line is observed in transcon- 
ductance and corresponds to a local maximum in the 
detector current /qpc (see also Figs. [7^ andjTjD). 

In order to explore the origin of this nonequilibrium 
feature, we will now discuss the pulse sequence step by 
step. The sketches in Fig. [5^ indicate the chemical po- 
tentials of the two QDs and its leads. In Figs. |4^ - 14]3 
the ground-state configuration at the black-white lines is 
0,1 with one electron charging QD C. This situation is 
indicated in sketch i of Fig. [5^. The short rise time of our 
pulses guarantee nonadiabatic transitions. Hence, at the 
very beginning of the square pulse, the electron is still lo- 
calized in the same QD as right before the voltage change. 
If pulsed onto the reconfiguration line, the configurations 
0,1 and 1,0 are degenerate, and thus, during the pulse, 
the electron tunnels with an interdot tunneling rate Fi 
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between the two QDs (sketch ii in Fig.|5p). Immediately 
after the pulse, the chemical potentials again correspond 
to the situation in sketch i. However, the new charge 
configuration depends on the pulse duration tp. Here we 
consider the incoherent limit where tp is long compared 
to the charge coherence time. In this limit and when 
^pTi ^ 1, we expect the two configurations 0,1 and 1,0 
to occur with equal probability (sketches iii-a and iii-b). 
However, for tpFi < 1, interdot tunneling is unlikely and 
the probability that the double QD stays in its ground- 
state configuration 0,1 is enhanced. 

In the case that the double QD occupies the excited 
configuration 1,0, it tends to relax towards its ground- 
state configuration 0,1 during the waiting time t^. Two 
possible relaxation channels 1,0 0,1 are indicated in 
sketches iv-a and iv-b in Fig. |5^. The first one includes 
interdot tunneling combined with energy relaxation of 
the electron, e. g. phonon-assisted tunneling (sketch iv-a). 
The other channel involves two resonant QD-lead tunnel- 
ing processes 1,0 0,0 with rate Fl followed by 0,0 
0,1 with rate Fr (sketch iv-b). For simplicity we summa- 
rize both relaxation channels 1,0 0,1 in the combined 
charge transfer rate F2 between two pulses — compared 
to the charge transfer rate Fi during a pulse. 

Charge transfer between QDs usually requires energy 
or charge exchange with the environment. In our case 
this results in a relatively small transfer rate F2. The 
transfer rate is much larger in the case of rapid resonant 
tunneling between the two QDs, which occurs if the con- 
figurations 0,1 and 1,0 are degenerate. Hence, pulses into 
this regime, namely onto the charge reconfiguration line 
or its linear extensions, yield Fi ^ F2. 

Fig. pp sketches the evolution of the probability P{t) 
for the (impulsed) ground-state configuration to be unoc- 
cupied in steady state, i. e. after many iterations of pulses. 
Our charge detector measures the time-average P of this 
signal. The combination of tp <C t^ and Fi F2 to- 
gether with F^^ > tw results in a nonequilibrium mean 
occupation appearing as a black- white line in Figs. - 

St. 



VI. QUALITATIVE DISCUSSION OF 
DIFFERENT COUPLING REGIMES 

In the case of Figs. - [4]:3 the pulse-induced black- 
white lines have exactly the same length as the charge 
reconfiguration lines. Here a non-equilibrium configura- 
tion apparently requires pulses directly onto the charge 
reconfiguration line (and not its extensions). This can be 
understood in terms of very asymmetric QD-lead tunnel- 
ing rates, namely a much larger tunneling rate between 
QD B and the left lead (Fl) compared to QD C and the 
right lead Fr <C Fl- During a pulse, from 0,1 onto the ex- 
tension of the reconfiguration line in region 1,1, a second 
electron might enter the double QD from the left lead, 
if Fl^p > 1. After the pulse, the transition 1,1 0,1 
brings the double QD back to its ground-state configura- 



tion on a time scale F^ < tp <C t^. Hence, the double 
QD spends most of its time in the configuration 0,1 which 
explains why the pulse-induced black-white line does not 
continue beyond the length of the charge reconfiguration 
line. On the other hand, during a pulse directly onto the 
reconfiguration line no charge is exchanged with the leads 
(sketch ii in Fig. |5^). If, after the pulse, the configura- 
tion 1,0 is occupied, the electron will very quickly escape 
from QD B to the left lead (sketch iv-b in Fig.jS^i). How- 
ever, refilling of QD C is much slower as Fr <C Fl and 
tunneling of an electron from the left lead into QD C is 
blocked. If, in addition, Fr^w < 1 we expect to find that 
the nonequilibrium configuration 0,0 is mostly occupied 
within the black-white line. This is indeed the case (as 
will be discussed in section VIII). 



Figures [6^ and [gJd show transconductance measure- 
ments for the same double QD as above — with very 
similar periodic pulse trains applied {tp = 0.5 ns and 
tw = 10 ns for Fig. [6|i) — but a set of slightly differ- 
ent gate voltages. Now the tunnel barrier between QD 
B and the left lead is tuned to be large compared to the 
barrier between QD C and the right lead Fl <C Fr (op- 
posite tendency compared to the data in Fig. [4|. In Fig. 

the pulse-induced black-white line is now extended by 
a dark line which merges into the (also dark) charging 
line 0,1 ^ 1,1 of QD B. Above this new line the charg- 
ing line of QD B is weakened. The data suggest that 
within the triangle, marked in Fig. [6^ by red corners, 
mostly the configuration 1,1 is occupied even though the 
ground-state there corresponds to 0,1. The sketch below 
Fig. [6^ explains this observation. The pulses starting 
from within the triangle bring the double QD into the 
region of the stability diagram with ground-state con- 
figuration 1,1. Even for Fl^,F^^ <C tp eventually the 
transition 0,1 1,1 occurs. The configuration 1,1 is not 
the ground-state between the pulses but still quite stable 
because the left tunnel barrier is large and the right QD 
is in Coulomb blockade. The observation of the discussed 
feature suggests Fr^p > Fl^w 

For Fig. [gJd the waiting time t^ = 20 ns is doubled com- 
pared to Fig. [6^ and the interdot coupling is increased, 
while all other parameters are left unchanged. On the 
one hand, the black-white line completely disappeared in 
Fig. [gJd, which suggests F2 ^ t~^ so that the ground- 
state is reoccupied quickly between two pulses. On the 
other hand, the just explained triangle remains present, 
because once the double QD is charged by two electrons 
the nonequilibrium state is again rather stable — despite 
of the increased interdot coupling. 

In Fig. |6fc {tp = 0.5 ns, t^ = 25 ns) the parameters are 
similar as for Figs. [4^ except that both tunnel barri- 
ers to the leads are now almost closed (Fl ~ Fr <C t~^). 
In this regime, charge exchange with the leads is absent 
and only the configurations 0,1 and 1,0 contribute to the 
nonequilibrium signal. Now, the black-white line is no 
longer a copy of the charge reconfiguration line but ex- 
tends throughout the region of ground-state configura- 
tion 0,1. 
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FIG. 6: Charge stability diagrams of a double QD for various 
QD-lead and interdot tunnel couplings (device and measure- 
ments as for Figs. |2] and [4). Sketches: relevant charge tran- 
sition processes (arrows) during pulses and between pulses. 
Crossed out arrows: energetically favored transitions which 
are very slow and, hence, result in nonequilibrium occupa- 
tions. Vertical lines depict tunnel barriers (smaller width 
equals larger tunnel coupling), (a-b) Asymmetric QD-lead 
tunnel couplings as indicated in the corresponding sketches. 
In a) tp = 0.5 ns and tw = 10 ns. In b) comparable param- 
eters except larger interdot tunnel coupling and tw = 20 ns. 
Within the triangle marked in a) the transition 0,1 1,1 
might happen during pulses, while the opposite transition 1,1 
0,1 is unlikely to happen between pulses, (c) Weak QD- 
lead tunnel couplings on both sides suppress the transition 
0,1 1,1 (tp = 0.5ns and tw = 25ns). (d) Weak QD-lead 
tunnel couplings on both sides combined with a very strong 
interdot tunnel coupling (tp — 0.5 ns and tw — 20 ns). 



Figure |6]i shows the extreme case of almost closed bar- 
riers to the leads (as for Fig. ^) but, in addition, very 
strong interdot tunnel coupling. Here, no effects of the 
applied pulse train (tp = 0.5 ns and tw = 25 ns) are vis- 
ible due to an extreme broadening of the transition 1,0 
^ 0,1. 



VII. RATE EQUATION MODEL 

In the following we develop a quantitative rate equa- 
tion model for the steady-state solution of the pulse cycle 



shown in Fig. ^jp. As an example we focus on the regime 
where the tunnel couplings to both leads are weak. Fig- 
ure ^ shows the stability diagram for such a situation 
in transconductance d/qpc/dV^ while Fig. [7^ plots the 
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FIG. 7: (a) Charge stability diagram of the same double QD 
as in Figs. [4^ - [4]: but for different QD-lead and interdot 
tunnel couplings. Plotted is the dc current ^/qpc after sub- 
traction of a linear background. A pulse train with tp = 0.6 ns 
and tw = 25 ns is applied, (b) (5/qpc along the cut through 
the stability diagram indicated in Fig. [7|l as a white dotted 
line (red triangles). The (black) filled circles show the mea- 
surement without pulsing, /lo, /qi and /oo are the current 
values expected for the charge configurations 1,0 and 0,1 and 
0,0 respectively, /peak marks the pulse-induced maximum of 
the detector current, (c) Relative detector current A// A/max 
(where A/ = /peak - /qi and A/max = /qo - /qi) as a func- 
tion of tp at tw = 25 ns for three different values of the gate 
voltage Vt2 (as indicated). Solid lines are model curves from 
Eq. [2] (see main text) (d) Filled circles are measured values of 
A//A/max as a function of tp and tw- The meshed surface is 
a model curve from Eq.|2]for the fit parameters Fj"^ = 0.74 ns 
and = 12.9 ns. Here we assume A// A/max = P (details 
in main text). 



dc detector current (5/qpc corrected by a linear offset. 
Instead of the black-white line observed in transconduc- 
tance, the absolute current maximum is now seen as a 
bright line. 

For simplicity, we consider purely exponential relax- 
ation processes. As above, interdot transitions 1,0 ^ 0,1 
are described by the rates Fi during pulses and F2 be- 
tween pulses. The validity of this approximation needs 
to be verified later, if, for example, the intermediate con- 
figuration 0,0 is involved (e.g. as for Figs. [4^ "[f]^)- To 
allow a direct comparison with experimental data we de- 
termine the probability, P(t), that the unpulsed ground- 
state configuration 0,1 is not occupied (either 1,0 or 0,0 
is occupied). Immediately before a pulse is applied (at 
t=0) and at the end of a pulse (at t = tp) P(t) can be 
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expressed by 



P(0) = P(tp)e-*-r. , 

P{tp) = 1 - [1 - P(0)] e-*pri . 



(1) 



After many repetitions, a steady state is reached, and 
P(0) = P{tp + tw) as indicated in Fig. Solving for 
P(0) and P(tp) using Eq|l]) leads to 



P(0) 



l-e*P^i 
l_e*wr2+tpri , 



Finally, integration over one pulse cycle in the steady- 
state limit results in the average probability that the un- 
pulsed ground-state is unoccupied 



/(I - [1 - P(0)] e-'^^)dt + yP{tp)e-'^^dt 





r2 



[P(0)-P(tp)]) . 



(2) 



VIII. QUANTITATIVE ANALYSIS AND 
SAMPLE CHARACTERIZATION 

Figure [TJ) plots SIqpq measured along the white dot- 
ted line in Fig. [7^ for two cases, namely, with (red tri- 
angles) or without pulsing (black circles). The average 
values of /qpc expected for the configurations 1,0 (/lo), 
0,1 (/oi) and 0,0 (/qo) are indicated. At the left end of 
the plot, /qpc — ho, which indicates occupation of the 
ground-state configuration 1,0. Next /qpc decreases to 
/o,i in two steps as the (split) charge reconfiguration line 
is crossed (see Figs. [4^ - c). For the case without puls- 
ing, /qpc stays approximately constant at /o,i before it 
increases to /o,o across the charging line between 0,1 and 
0,0. 

In contrast, the curve produced while pulsing contains 
an additional maximum labeled /peak within the region 
of the ground-state configuration 0,1. This maximum 
corresponds to the bright line in Fig. [7^. Its amplitude 
provides not only information about the average proba- 
bility P, but also about the relevance of the two possible 
relaxation channels depicted in the sketches iv-a and iv-b 
in Fig. [5^. For the data presented here, we observe that 
/peak can exceed /i^o and even takes values up to /o,o- 
This behavior contradicts the expectation for a combi- 
nation of incoherent resonant tunneling (sketch ii in Fig. 
Isk) and direct energy relaxation via interdot tunneling 
(sketch iv-a), because the overall occupation of one elec- 
tron in the double QD is preserved in these processes. 
In the case of incoherent interdot tunneling, /peak would 
be limited by 0.5 x /iq. Hence, our system favors the 
other relaxation channel 1,0 ^ 0,0 ^0,1 (sketch iv-b in 



A// A/max, where A/ 
/oi, as a function of tp for fixed t 



Fig. |5^) which involves charge exchange with the leads. 
Moreover, /peak ^ /i,o implies that in this case the in- 
termediate configuration 0,0 is mostly occupied, which 
requires that the tunneling process 1,0 ^ 0,0 is faster 
than the subsequent transition to the ground-state 0,0 
^0,1. In this limit the two-step tunneling process can 
be approximated by a single exponential decay rate T2 
and our model expressed in Eq. Q can indeed be ap- 
plied. However, it is possible to tune the double QD into 
different regimes (compare Fig. [6|, some of which would 
require modification of the rate equation model. 

Figure ^ presents the relative current change 

/peak — /oi and A/max = /qO " 

25 ns. The three dif- 
ferent symbols correspond to three different voltages ap- 
plied to gate t2 which controls the interdot tunneling rate 
Fi. As expected A//A/max grows faster for a smaller in- 
terdot tunnel barrier t2 (less negative Vt2) and a larger 
tunneling rate Fi. As long as the occupation time of con- 
figuration 1,0 is very short compared to that of 0,0 our 
model allows the direct interpretation A//A/max = P 
and, therefore, a quantitative comparison with Eq. ([2|. 
The lines are the according model curves for F^^ = 75 ns 
and = 0.39 , 0.92 and 2.33 ns for Vt2 = -285 , -295 
and -300 mV. 

A more complete comparison between similar data 
and our rate equation model is presented in Fig. [7]i, a 
plot of P as a function of tp and t^. The measured 
data (closed circles) are overlayed with a two-dimensional 
model curve representing Eq. (|2| with the fit parameters 



0.74 ns, r:; 



12.9 ns 



The overall smooth 



increase in P with increasing tp and decreasing indi- 
cates incoherent tunneling between the two QDs in these 
measurements. 



IX. SUMMARY 

In conclusion we have demonstrated how charge spec- 
troscopy is capable of resolving single electron charging 
events in complex nanostructures at time scales below 
1 ns. These experiments require a large bandwidth and 
it is crucial that the sample holder and other rf com- 
ponents are well designed to avoid refiections and cross- 
talk between control gates. In our measurements we have 
combined periodic trains of ultra-short voltage pulses on 
specific gates with a continuous (dc) charge detection us- 
ing a QPC. A rate equation model can be used to quanti- 
tatively analyze the data and to separate tunneling rates 
from energy relaxation times. Different relaxation chan- 
nels in the system can be identified and controlled. This 
pulsed-gate technique, demonstrated here for a full char- 
acterization of a double QD, can be used for coherent 
control of charge qubits or coherent transport by adia- 
batic passage [29l |30] . Our approach promises scalability 
towards an array of qubits in which each qubit is sepa- 
rately addressable while mutual qubit entanglement can 
be fully controlled. 
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